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We investigate cyclotron motion in graphene monolayers considering both the lull quantum dy- 
namics and its semiclassical limit reached at high carrier energies. Effects ol zitterbewegung due to 
the two dispersion branches ol the spectrum dominate the irregular quantum motion at low energies 
and are obtained as a systematic correction to the semiclassical case. Recent experiments are shown 
to operate in the semiclassical regime. 
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I. INTRODUCTION 



our results with the Shubnikov-de Haas expriments by 
Novoselov et alA. Section HVl contains our conclusions. 



The experimental discovery of graphene. i.e. single lay- 
ers of graphite, has rapidly led to an extraordinary vast 
and still growing interest in this material, both exper- 
imental and theoreticalii^'^. Compared to conventional 
two-dimensional electronic systems, the peculiar proper- 
ties of graphene mainly stem from its linear dispersion 
near the Fermi energy, and the chiral nature of electronic 
states entangling the momentum and sublattice degree 
of freedom^ii^. In particular, studies of carrier transport 
in perpendicular magnetic fields have revealed unusual 
features like a cyclotron mass being proportional to the 
square root of the particle density^, and, most spectacu- 
lar, a quantum Hall effect occuring at half-integer filling 
factors^i^. Another partially related feature of graphene 
is its spectrum of Landau levels being non-equidistant in 
energy and proportional to the square root of the mag- 
netic field,— i^i^iiiSii^, properties also not shared by usual 
two-dimensional electron gases in semiconductor struc- 
tures. 

In this communication we analyze the cyclotron mo- 
tion in graphene considering both the full quantum 
dynamics as well as its semiclassical limit. The lat- 
ter findings can directly be compared with the the re- 
sults of Shubnikov-de Haas measurements reported on in 
Refi^. We also comment on a recent preprint by Rusin 
and Zawadzki studying cyclotron motion in graphene 
which appeared while the present work was reaching 
completioni^. These authors concentrate on the full 
quantum dynamics using a numerical approach anal- 
ogous to very recent workii on two-dimensional elec- 
tron gases in semiconductors, and they stress the role 
of zztterfcewe^Mn^iHai^aeiii^as^^^^^^. As we 
shall see below, signatures of zitterbewegung can also be 
seen in semiclassical corrections to the classical limit of 
the underlying quantum dynamics. 

In the following section |TT] we investigate the full quan- 
tum dynamics of the system using both a numerical 
and an anyltical approach. The latter one is based on 
the analogy of the Hamiltonian to the Jaynes-Cummings 
model steming from quantum optics. This analogy is the 
starting point for the analysis of the semiclassical limit 
to be discussed in section IIIIl There we also compare 



II. QUANTUM DYNAMICS 

For a graphene sheet in a perpendicular magnetic field, 
the single-particle states around one of one of the two in- 
equivalent corners of the first Brillouin zone are described 



(1) 



with (using standard notation) tt = p + eA/c and v ~ 
lO^m/s. The Pauli matrices describe the sublattice or 
pseudospin degree of freedom, and the Zeeman coupling 
to the physical electron spin has been neglected. The 
label = ±1 determines which corner of the Brillouin 
zone in considered; in what folows we shall concentrate 
on r" = 1 The Heisenberg equations of motion read 



-JI'^hW = -p<yH[t) X e, 
— crif(t) = -^T^H{t) X anKt) 



(2) 
(3) 



where t ~ ^ hc/\eB\, {—e)B > 0, is the the magnetic 
length, and is the unit vector along the z-direction. 
The magnetic field B = V x A = Bcz is assumed to 
point along the negative z-direction, B < 0. The position 
operator f = (x, y) can be given in terms of the kinetic 
momentrum tt via the usual relations 



y = yo 



eB 



(4) 
(5) 



where ro = {xo,ya) is conserved, [?i,ro] = 0. 
For cyclotron motion in a usual non-interacting two- 
dimensional electron gas, the vector rg describes the cen- 
ter of the classical circular orbits. As we shall see below, 
this is also the case for the classical limit of cyclotron 
motion in graphene. 

Defining the usual bosonic operators 
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V2h 



(tTj; 4- ITTy) 



(a)^ 



(6) 
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fulfilling [a, a+J = 1 the Hamiltonian reads 



hv 



n='-^V2{aa- +a+a+) 



(7) 



where cr± = (cr^±iCT2')/2. The energy scale of this Hamil- 
tonian is given as a function of the magnetic field by 



— w 2omeV 



Tesla 



(8) 



while its length s cale is the usual magnetic length, £ = 
257A/^/B/^es[a.. The well-known eigenstatea^ of the 
Hamiltonian ([T]) are given by |0,t) with energy Eq ^ 
and. for n > 0, 



|n,±) = ^(|n,T>±|n-l,i» 



(9) 



with energy = ±{hv / £)\/2ri. Here n is the Landau 
level index, and the arrows are obvious standard notation 
for the sublattice spin states. 



A. Numerical Approach 

The operator- valued equations of motion ([2]), ([3]) do 
not seem to allow for a full exlicit solution. However, it 
is straightforward though somewhat tedious to numeri- 
cally evaluate the time evolution of momentum, position, 
and spin operators. A similar approach was performed 
recently in Ref.— investigating cyclotron motion in semi- 
conductor quantum wells with spin-orbit coupling. In 
fact, the present case of graphene is technically clearly 
simpler than the previous one and was also studied very 
recently in Ref.—. For such numerical simulation it is 
convenient to work in the Landau gauge A = (0, Bx, 0) 
with the initial state lip) being a direct product of an 
orbital and a spin state. 



(10) 



where the spinor components are related to the usual 
polar angles (p of the initial pseudospin direction via 
K = exp(— z(p/2) cos(i9/2), A = exp(z</j/2) sin(i?/2). As a 
generic initial orbital state we consider 



r+ikoy 



(11) 



i.e. a normalized Gaussian wave packet of spatial width d 
and initial momentum fifco along the y-axis, i.e. the direc- 
tion of translational invariance of the Hamiltonian. The 
initial position of the particle is at the origin, {tp\f\tlj) — 0, 
and its energy is given by 

E ^ {ij\n\tp) ^ hvko^ {kX - kX) (12) 



with a quantum mechanical uncertainty of 



{Any 



hv 



PI ( rf2 + 2^2 (I 



+kl£'^ {l+{RX^nXf^ j (13) 



Note that the initial state ([TU)) has in general non- 
vanishing overlap with single-particle eigenstates of the 
form ^ of both positive and negative energy. As to 
be discussed below, this fact leads to additional oscil- 
lations in the time evolution that can be viewed as 

We emphasize that, by the very construction of the 
Hamiltonian ([T]) , the wave number fco is to be interpreted 
relatively to the wave vector of the chosen corner of the 
first Brillouin zone^. The authors of RefJ^, however, 
consider a time evolution of an initial state of the form 
pop under the simultaneous action of the Hamiltonians 
of both inequivalent corners of the Brillouin zone, a theo- 
retical modeling whose physical meaning remains rather 
unclear. Moreover, these two Hamiltonians are assumed 
to differ just in a global sign, i.e. one Hamiltonian is the 
negative of the other, which is at odds with the micro- 
scopic tight-binding description of graphene^. 

For an infinite sheet of graphene it is natural to con- 
sider initial conditions where both sublattices have the 
same quantum mechanical weight, |k| = |A|, i.e. the sub- 
lattice or pseudospin lies initially in its a;?/-plane. Fig. [1] 
shows several examples for numerically evaluated trajec- 
tories {rnit)) ■= (V'I^-h(OIV') with this type of initial 
condition. For all further details of these conceptually 
straightforward but technically somewhat involved nu- 
merical simulations we refer to RefsJ^iii. In the two top 
panels of Fig. [T] the sublattice spin is initially coUinear 
with the momentum. These simulations can be com- 
pared with Fig. 5 of Refii^ where the authors consider 
a time evolution under a single Hamiltonian (not two of 
them) given by Eq. ([T]). Indeed, all results given in that 
figure are essentially reproduced by our own simulations 
underlying the present work. 

The remaining panels of Fig. [1] contain simulations 
where the sublattice spin is in the initial state not 
coUinear with the momentum. Note that the initial ve- 
locity (vh) = (th) is determined by the initial direction 
of the sublattice spin via 



(14) 



Moreover, further simulations show that there is no dy- 
namics in the initial direction of the momentum, i.e. 
iUHit)) = 0, if the sublattice spin is initially perpen- 
dicuUar to it, i.e. in the xz-plane. In such a case also 
no sublattice spin component coUinear with the initial 
momentum develops in the time evolution, (cr^(t)) = 0. 
These observations were partially already made in Refii^ 
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and can be understood from the equations of motion ([2]), 
([3]). In general, the trajectories seen in Fig. [T] are rather 
irregular which is a consequence of the many different 
excitation frequencies being present in the spectrum at 
low energies. In RefJ^ contributions to the time evo- 
lution involving transition frqucncies between states of 
positive and negative energy have been regarded as an 
effect of zitterbewegung. Indeed, as to be shown below, 
such type of zitterbewegung also occurs as a correction to 
the classical limit. 

As already mentioned in RefsJ^iii the numerical ap- 
proach discussed above is technically limited to initial 
states having significant overlap with rather low Landau 
levels only, i.e. the method is restricted to the regime 
dominated by quantum effects. In the following we shall 
therefore explore the semiclassical limit using an analyt- 
ical approach. 



B. Analogy to the Jaynes-Cummings model 

Further analytical progress regarding the full quantum 
dynamics can be made by exploiting the fact that the 
Hamiltonian ([7]) is formally equivalent to the Jaynes- 
Cummings model for atomic transitions in a radiation 
field. A similar observation has been made recently in 
Refi^ for the two-dimensional electron gas in semicon- 
ductor quantum wells with Rashba spin-orbit coupling. 
The Jaynes-Cummings model has been studied very in- 
tensively in theoretical quantum optics, and the time evo- 
lution of orbital and spin operators has been obtained in 
terms of analytical but rather implicit expressions2Si21. 
Using the method described in Refsj^i^l one can solve 
for the time-dependent position operators in the Heisen- 
berg picture as 



XH{t) + iyH{t) = xo + iyo 



Ul- — L0+ 



(uj^ {Kx + IKy) i - S^O'^) 
(w+(K:,-HiK.y)^-2^Cr+) , (15) 



with K = tt/S, and the operator-valued frequencies lo±_ 
are given by 



-n±\ + 2 



(16) 



All operators on the r.h.s of Eq. (|15p are in the 
Schrodingcr ppicture, i.e. at time t = 0. As in the case of 
the two-dimensional electron gas with Rashba spin-orbit 
coupling investigated prcviouslyii*^, the result (fT5|l is 
still rather implicit and difficult to evaluate for a given 
initial state, mainly due to the operator character of the 
frequencies lu±. However, as we shall see in next sec- 
tion, the above result provides a very natural access to 
the classical limit of the dynamics including semiclassical 
corrections. 



III. SEMICLASSICAL LIMIT 

Cyclotron motion of massful electrons is in the simplest 
case just described ba a kinetic (effective-) mass term, and 
the quantum mechanical and the classical equations of 
motion coincide, the latter ones being the obvious and 
well-defined classical limit of the former. However, the 
classical limit of the Hamiltonian ([T]) describing mass- 
less fermions appears prima vista not as obvious. In 
general, the classical limit of a quantum system is ap- 
proached in the limit of high energies. For our problem 
here this means that the energy E = (H) of electron 
must be large compared to the charcteristic energy scale 
hv/£ of the Hamiltonian, which is equivalent to the con- 
dition (k • (t)^ 3> 1. In what follows we shall assume an 
initial state of the form (jlOp with momentum and sublat- 
tice spin being parallel to each other. This leads to the 
condition 



ko£ > 1 



(17) 



This very natural classical limit of the graphene model 
in a perpendicular magnetic field corresponds to the 
"strong-coupling scenario" discussed in Refi^. We note 
that a similar limit is reached for negative energies of 
large modulus. Here momentum and sublattice spin are 
initially antiparallel. 

In the limit of large energies E, the operator- valued 
frequencies w± can be replaced with classical variables, 



-E± \ E^ 



(18) 



Then, treating also the the quantites r, fo as well as k, a 
as classical variables (not as operators) one derives from 
the full quantum result in the above classical limit 



x{t) - Xo 



£ 



LU- (Kx£sm{iU+t) — Ky£ cos{Lu+t)) 



+ - {—a^ sin{uj+t) + cos((jj+t)) 
— W-i- {Kx£sm{uj-t) — Ky£cos{u!-t)) 

V 

— — {—(T^ sin(w_t) + cr*' cos(cj_f)) 



(19) 



yit) - yQ 



£ 



LU- {kx£ cos{LLi+t) + Ky£sm{Lu^t)) 



-- {a^ cos([j+<) -I- 17^ sin(cj_|_i)) 
-a;+ [nxi cos{ll! -t) + Ky£sin{u!-t)) 



+ - {(7^ coa{uj_t) + sin(a;_t)) 



(20) 
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Again the quantities k, a in the rectangular brackets 
on the r.h.s. are at time i = 0. Moreover, in the 
classical limit E ^ fmjl we have hw^ w —2E and 
fia;+ « (hvjti^jE. It is instructive to rewrite the lat- 
ter expression in the usual form of a cyclotron frequency 
hwj^ Wc = \eB\lmcC. Here the cyclotron mass 
is given by the well-known scmiclassical expression^i^^ 
771c = {h^ /2i:)dS/dE ~ Ejv^ where S = Trfcg is the area 
enclosed by a cyclotron orbit, and E = hvko. Now, ex- 
panding the above expressions in first order in (hv/£)/E 
one finds 

x{t) — Xq — Kxi^ Sm{uJct) — Ky£'^ cos{uJct) 

+ ^ (—f^^ sm{ujct) + (7^ cos(a;cO 
E V 

-a"" sm{2Et/h) - cos{2Et/h)^ , (21) 
y{t) — Ho = Kxi^ cos{u}ct) + Ky£^ sin(a;ct) 




(T^ cos(ci;ci) + o-^ sin(cjct) 



-cr^^ cos{2Et/h) + sin{2Et/h)j . (22) 

The first lines of the above r.h.s. describe the classi- 
cal cyclotron motio n^^i^^ with an energy-dependent cy- 
cloron mass^. The other contribution are lowest-order 
semiclassical corrections to this classical limit. In partic- 
ular, the terms in the last lines oscillate with the large 
frequency 2E/h ^ ujc which equals the energy separa- 
tion 2E = 2hvko between states of positive and negative 
energy at given wave vector in the absence of a mag- 
netic field^. Therefore, these semiclassical correction can 
be viewed as an effect of zitterbewegun g^^^^^^^^ . We note 
that in connection with graphene thne term "zitterbewe- 
gung " was also used recently to express the fact that the 
velocity operator (fT4|) fails to commute with the Hamil- 
tonian (P)^. 

So far we have concentrated on semiclassical dynamics 
at large positive energies E ^ hv/l. In the analogous 
case of negative energies of large modulus one obtains a 
similar result with the frequencies uj± being interchanged, 
cij+ « 2|_B|, Lj„ « — Wc. 

Finally we note the equivalence of the following three 
conditions: 

E->hiUc = (23) 
^ E^^ ihvko)' » (f )' (24) 
^ (M)^>1 ■ (25) 

Thus, the condition ((23|) is fulfilled whenever (fT7|) is valid. 



On the other hand, the condition ((T7|) can be rewritten 
as Tc ^ ^ with Tc = kgi'^ being the classical cyclotron ra- 
dius. This is indeed the usual textbook criterion for the 
validity of semiclassical approximations to cyclotron dy- 
namics in solids^^ and confirms again our above strategy 
of obtaining the semiclassical limit. 

Let us now compare our results with the Shubnikov- 
de Haas measurements by Novosclov et alA. The data 
presented in Fig. 2a of Ref.— was obtained at an elec- 
tron density of ti sa 4.3 • lO-'^^cm^ and magnetic fields of 
up to about B ~ lOT, corresponding to a Fermi wave 
vector of fcj = a/ttu w 3.7 • lO^cm"^ and a magnetic 
length of £ > 8lA. Thus, kfl > 3, and the criterion 
(fT7)) is still reasonably fulfilled even for the highest fields 
used in those measurements. Correspondingly, we have 
{hv/t)/E < 0.33 with E = hvkf, showing that the semi- 
classical correction obtained in Eqs. pT|) . ((22|) are suf- 
ficiently suppressed which allows for an interpretation 
of the experimental data in purely classical terms. At 
higher magnetic fields, quantum effects dominate, and 
quantized Hall transport is observed^^^. 



IV. CONCLUSIONS 

We have investigated cyclotron motion in graphene 
monolayers considering both the full quantum dynam- 
ics and its semiclassical limit obtained for high carrier 
energies. At low energies, the quantum dynamics leads 
to rather irregular particle trajectories dominated by ef- 
fects of zitterbewegung due to the two dispersion branches 
of the spectrum. The semiclassical limit of the system is 
obtained using the analogy of the Hamiltonian with the 
Jaynes-Cummings modeliii^. A similar analysis of the 
classical limit can be done for cyclotron motion in semi- 
conductor quantum wells with spin-orbit interactionii*2^. 
The semiclassical limit of cyclotron dynamics in graphene 
is described by the usual cyclotron frequency and the 
characteristic frequency of zitterbewegung which occurs as 
a semiclassical correction. Recent experiments are shown 
to operate in the semiclassical regime. 
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FIG. 1: Orbital dynamics of a wave packet of initial width 
d — 1.0£ and group wave number ko = 1.0/£. In all cases 
the sublattice spin lies initially in the ay-plane, i.e. — n/2 
leading to \k\ = \X\, and the total simulation time is always 
50£/v. In the left and right top panel, the sublattice spin 
is initially coUinear with the momentum with tp = n/2 and 
= 3n/2, respectively. In the middle panels we have <^ = 7r/4 
(left) and if = Bn/A (right) as initial conditions, while in the 
bottom panels ip = Sn/A and (p — 7tt/4: were used 



